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Abstract: The enumeration of BPS bound states in string theory needs refinement. 
Studying partition functions of particles made from D-branes wrapped on algebraic Calabi- 
Yau 3-folds, and classifying states using split attractor flow trees, we extend the method 
for computing a refined BPS index, [1]. For certain D-particles, a finite number of mi- 
crostates, namely polar states, exclusively realized as bound states, determine an entire 
partition function (elliptic genus). This underlines their crucial importance: one might 
call them the 'chromosomes' of a D-particle or a black hole. As polar states also can be 
affected by our refinement, previous predictions on elliptic genera are modified. This can 
be metaphorically interpreted as 'crossing-over in the meiosis of a D-particle'. Our results 
improve on [2], provide non-trivial evidence for a strong split attractor flow tree conjecture, 
and thus suggest that we indeed exhaust the BPS spectrum. In the D-brane description of 
a bound state, the necessity for refinement results from the fact that tachyonic strings split 
up constituent states into 'generic' and 'special' states. These are enumerated separately 
by topological invariants, which turn out to be partitions of Donaldson-Thomas invari- 
ants. As modular predictions provide a check on many of our results, we have compelling 
evidence that our computations are correct. 

Keywords: Elliptic genus, Calabi Yau, Topological String, Donaldson-Thomas 
invariant. Wall-crossing. 



Contents 



1. Introduction 1 

2. Exact enumeration of BPS microstates of a D particle 5 

2.1 An index for BPS bound states using Donaldson-Tliomas (DT) invariants 6 

2.2 Polar states: tlie cliromosomes of D-particles / black iioles 7 

3. Elliptic genera from split flow trees and DT invariants 9 

3.1 Polar states on the sextic liypersurface in VKP^;^]^^2 9 

3.2 Polar states on the octic hypersurface in W^P^-^j^^^ 11 

4. Reflned enumeration of bound states and elliptic genera revisited 14 

4.1 The main idea: non-trivially fibered moduli spaces of bound states 14 

4.2 Algebraic techniques to deal with special constituent states 16 

4.3 Special points and DT partitions J\fj^r^^\o, 2) 17 

4.4 Special curves and DT partitions 7V^?^''^(1, 2) 22 

5. Reflned predictions for elliptic genera 24 

6. Discussion 31 

6.1 Absence of scaling solutions? 31 

6.2 The meiosis of D-particles/black holes 33 

A. Setup and conventions 36 

B. Modifled elliptic genera and D particles 37 

C. The split attractor flow tree conjecture 39 

D. Methods to establish split flows outside the supergravity regime 42 



1. Introduction 

By virtue of the correspondence principle, BPS black holes in 4d and 5d string theory 
compactifications can be studied as bound states of D~brane systems. These D~branes 
are wrapped such that they form a particle from a 4d (or 5d) point of view: a so-called 
D-particle. If the D-branes carry a sufficient amount of charge, the system will form a 
black hole through the gravitational backreaction as soon as the string coupling is turned 
on. Also a D-particle carrying minimal charge is a very interesting toy model, which offers 
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itself for studying the basic mathematical principles governing BPS black holes, although 
the direct interpretation as a black hole is lost. 

In this paper, the partition functions of D~particles modeled as mixed (magnetic/electric 
charge) ensembles with branes wrapped on algebraic Calabi-Yau 3-folds are studied. The 
microstates of key interest are so-called polar states, which are given the interpretation of 
chromosomes of D-particles or black holes in subsection 2.2. One of the main virtues of 
polar states of black holes is that they are always multi-centered solutions. In the case of 
a D-particle, polar and some non-polar bound states are enumerated, and we propose a 
refined computational scheme to calculate indices for these bound states. At the same time, 
we split up Donaldson-Thomas invariants, used to enumerate constituents of these bound 
states, into Donalds on- Thomas partitions. The index for which we propose a refinement 
was used in [3] . The authors of this paper also state that a moduli space of a bound state 
is in general a fibration. This was examined for an example in [1] and a refinement of the 
index enumerating a bound state was suggested there. This work can be seen as a further 
development of these latter ideas. 

Interest in mult i centered BPS black hole solutions 

In recent years, great interest has arisen in multi-centered black hole solutions within 
the BPS spectrum of type II string theory compactifications. This has several reasons, of 
which we shall list a few of special interest for the development of our work: 

• According to the BPS black hole attractor mechanism, a spherically symmetric 
(single-centered) black hole solution can be found in A/" = 2, d = 4 supergravity 
theories, independently of the chosen values of the vector multiplet scalar fields at 
infinite distance of the black hole too := t{r = oo), where t{r) denotes the value of the 
scalars in dependence of a radial spacetime coordinate r. The value of these moduli 
will be driven to the so-called 'attractor value' at the event horizon, t^, := t{r = 0). 
The values of these moduli form a line in moduli space (connecting too and t^, ) , which 
will be referred to as a single flow. G. Moore put forth a correspondence between 
spherically symmetric solutions to the BPS attractor equations and BPS states in 
string theory. Studies by G. Moore in 1998 [4, 5] showed that this correspondence 
does not hold. For a full correspondence, (at least) stationary multi-centered solu- 
tions must also be considered. 

• The connection between the supergravity description of BPS black hole solutions and 
the D-brane description has received completely new impulses from the studies on 
multi-centered solutions and bound states over the last decade. In [6] it was shown 
how a D-brane system wrapped around a compact manifold, yielding a single point 
in space, can transform into a bound state of two (or more) constituents, in the 
gravitational description. Genuine bound states of black holes, subject to specific 
equilibrium distances for stability, were discovered by Denef and collaborators [7, 8]. 
A review how to construct these multi-centered solutions was given in [9]. By now, it 
seems fair to say that it has been acknowledged that they form a very prominent part 
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of the spectrum. A sample study of the BPS spectrum for a type II compactification 
on the quintic 3-fold provides a nice illustration [10]. 

• The formulation of the OSV conjecture [11] released a period of very active research 
on the connection of black hole entropy and topological string theory; see [12] for a 
review. Basically, the conjecture suggests a relation of a black hole partition function 
with the topological string partition function of the form, Zbh ~ |-2^top|^- Attempts 
to prove this conjecture have proven to be extremely difficult. One attempt by Denef 
and Moore, which is of special relevance for the present work, put multi-centered 
black hole solutions to the center of the stage [3]. The authors concretized the 
conjecture using a D4-D2-D0 partition function as a black hole partition function, 
built from a mixed ensemble, keeping the D4-charge fixed and varying over the 
D2/D0 charges. This partition function displays modular invariance and it turns out 
that it is completely determined by knowing the degeneracy of a finite number of 
microstates: the polar states^. Such a relation was derived by the authors of [3] and 
independently in [13], subject to certain amendments [14]. These polar states are 
exclusively realized as multi-centered solutions, as will also be discussed in section 2. 
A large part of the work presented in this paper consists of the enumeration of polar 
states for D-particles. 

• Another finding of [3] is, that upon scaling up the charge of a brane system, the 
entropy of multi-centered microstates carrying that charge will always start to dom- 
inate over the entropy coming from single-centered states. This is puzzling, as the 
entropy of multi-centered solutions scales with the cube, 5 — )• A^5 when scaling D- 
brane charges Q — )• AQ, and not with the square, as is the case for single-centered 
solutions, S — )• A^S. For a (large) black hole, that carries large charges, this implies 
that one would expect a far too large entropy and this puzzling discovery has been 
given the name entropy enigma. This again stresses the importance of multi-centered 
solutions. A possible connection exists with another phenomenon discovered, namely 
the appearance of scaling solutions. Scaling solutions appear as multi-centered BPS 
solutions, which however allow the continous variation of a parameter, the distance 
between constituents. This implies that one can vary this 'scaling modulus' and 
choose two centers to lie so close together in spacetime that their throats virtually 
melt together, and the constituents become indistinguishable to an outside observer. 
The precise meaning of these solutions is at this moment still a point of ongoing 
discussion, and their implications for the split attractor flow conjecture (discussed 
in section C) are not completely clear. These issues have been recently addressed 
in [15, 16, 17]. We will discuss our results in the light of these issues in a separate 
section, 6.1. 

Just like single-centered BPS black holes, multi-centered BPS black hole solutions are 
also governed by an attractor mechanism. It was originally extended to multi-centered 

^These polar states are typically states with a low D2/D0-charge; upon increasing the D2/D0 charge, 
the states become non-polar at some point. 
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black holes in [8]. The image of the scalar moduli in moduli space forms a line, that splits 
(possibly several times) and runs from a background value ioo to a split point (and maybe 
further split points), in order to end at two (or several) attractor points, tu,t2*,---, one 
for each center of the bound state. This image has been given the name split attractor 
flow tree and has been conjectured to be an existence criterion for a multi-centered black 
hole solution. The meaning of split flow trees extends beyond supergravity, as has been 
confirmed by studies on the quiver description of D-brane systems in [6]. In fact, the split 
attractor flow tree conjecture states that (single and) split flows completely classify the 
BPS spectrum of type II string theory. The results in this paper can be interpreted as 
strong evidence for this conjecture. 

Main goal: index refinement for BPS bound states 

In [1], a study very similar to those in this paper was performed, and a refined prescription 
for calculating an index for a BPS bound state was proposed. The reason for the necessity 
of this refinement lay in the fact that the tachyon fields connecting the two constituents of 
a bound state did not perceive all the individual microstates of a constituent generically. 
The investigations in this paper go in the same direction, however the reason why the 
tachyon fields do not perceive all constituent states of a bound states generically, will be 
different. This leads us to new techniques to calculate the refined index. Non-trivial checks 
on our method yield exact confirmation of the predictions from modularity on degeneracies 
of certain states, suggesting that our technique is indeed correct, and beyond that, that 
the strong split attractor fiow conjecture might well be completely accurate. Whereas the 
'naive' index to enumerate a bound state is a simple product, we argue, as in [3, 1], that 
the moduli space is rather a fibration. This means that the index falls into a sum of several 
pieces, which can be grouped logically. This leads us to distinguish between those states 
of a constituent, in our case a D6-brane system enumerated by a Donaldson-Thomas in- 
variant, that are perceived generically, and those that are perceived as special states. This 
also allows the definition of Donaldson-Thomas partitions, which enumerate those generic 
and special states separately. The main idea of our technique will be explained at the 
beginning of chapter 4, which also contains our results on the refined bound state index 
and on Donaldson-Thomas partitions. 

Organization of this paper 

• Section 2 covers some basic material and may be skipped by the expert reader. Sub- 
section 2.1 briefly recapitulates the index developed for the enumeration of BPS 
bound states in type II string theory compactiflcations in [3], extensively used and 
eventially improved on in this paper. For the sake of stressing their importance, 
special bound states, namely polar states, which determine elliptic genera through 
modularity, are given the interpretation of chromosomes of a D-particle / a black hole 
in subsection 2.2. The reader unfamiliar with techniques used in this paper will flnd 
short, but more or less self-contained appendices covering the relevant background 
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material used in this paper. We mostly follow the setup of [3] and [1] . The reader will 
find more explanation on the setup of this paper, in particular D-branes wrapped 
on algebraic one-modulus Calabi-Yau varieties, in appendix A, a short account of 
(modified) elliptic genera in appendix B, and a short introduction to split attr actor 
flow trees as well as a formulation of the split attractor flow tree conjecture is pre- 
sented in appendix C. Eventually, to concretely establish the existence of BPS states 
for our study models, we will use numerical techniques, involving the use of mirror 
symmetry, adapted from [10]. This is necessary, as we work with very small charges, 
and therefore instanton corrections to the central charges of brane systems become 
dominant. More details on these techniques can be found in appendix D. 

• Section 3 presents our results for polar states on two Calabi-Yau's, described as 
hypersurfaces in weighted projective spaces. These results allow a prediction on the 
partition function for a D-particle, and exactly match the results in [2]. 

• Section 4 presents our most interesting findings. For some tractable non-polar states 
for two of our models, calculations based on the 'naive computation' of a BPS index 
show a discrepancy from the result predicted by modularity. We demonstrate that our 
computations based on a refined index for BPS bound states however show an exact 
matching with the prediction. At the same time, we use our procedure to distinguish 
between constituent states that are perceived by the tachyon field generically, and 
those that are perceived differently: we use the terms 'generic' and 'special' states 
and call the invariants that enumerate these D6-brane systems, Donalds on- Thomas 
partitions, J\f^^{(3,n),J\fj^^{l3,n), where the superscripts 'g' and 's' stand for 'generic' 
and 'special', respectively, and where A'dtI/?,?^) = A/'Q^(/3,n) + J\fj^r^{l3 , n) . 

• In section 5, we show that the refined prescription to compute BPS indices for bound 
states also alters the enumeration of polar states, leading to a new prediction for the 
elliptic genus for a degree ten hypersurface, embedded in a weighted projective space. 
This prediction is slightly different from the results in [2]. 

• Section 6 is a discussion of our results. We devote some time to the question whether 
anything can be learned about a 'tentative' part of the classical BPS spectrum of a 
D-particle, namely scaling solutions, in a separate subsection, 6.1. In another subsec- 
tion, 6.2, we interpret our most interesting results, the necessity for the refinement of 
bound state indices, as an artefact of a 'sort of meiosis for D-particles / black holes'. 
We conclude with some implications and possible directions for future research. 

2. Exact enumeration of BPS microstates of a D— particle 

We model D-particles using mixed ensembles of D4-D2-D0 branes wrapped on a Calabi- 
Yau 3-fold X, which we choose to be a hyperplane in a weighted projective space. A 
D4-brane with charge p = 1 is wrapped on the hyperplane class divisor, and kept fixed, 
while varying over all possible U{1) worldvolume fiuxes F G H'^{X) and various numbers 
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of bound DO-branes. We follow the setup of [3, 1]. The reader can find more details on 
our setup and conventions in appendix A. 

The existence of D4-D2-D0 BPS states will be inferred from the the existence of split 
attractor flow trees. The reader unfamiliar with this technique can find a short introduction 
in appendix C. As we inspect low-charged BPS states, higher order curvature corrections 
to the action become important, and central charges of brane systems will have to be 
calculated by exploiting mirror symmetry. Details on this procedure can be found in 
appendix D. 

2.1 An index for BPS bound states using Donaldson Thomas (DT) invariants 

Our goal is to describe D4-D2-D0 configurations as bound states of D6-D2-D0 and D6- 
D2-D0, possibly carrying U(l) fiuxes, whereby the D4 charge is induced by the latter 
fluxes. This will allow us to factorize the indices of D4 systems as products of D6 and D6 
indices. In order to compute the BPS indices of the D6 and D6 systems, we will make use 
of Donaldson-Thomas invariants. An invariant N-[)T{f3,n) computes the Witten index of 
a system with a D2 brane wrapping a curve of homology class /3, and a collection of DO's, 
such that the total DO charge equals n. Although the U(l) flux on the D6 interacts with 
these lower branes, it does not alter the Witten index. In mathematical terms, the DT 
invariants compute the dimensions of the moduli spaces of the ideal sheaves corresponding 
to curves and points on the Calabi-Yau. They are indeed conjectured [18, 19] to contain 
equivalent information as the Gopakumar-Vafa invariants [20, 21], which count the states 
of M2 branes with momentum, where the M2's are wrapped on holomorphic curves. By the 
conjectured identity between the generating functional for GV invariants and DT invariants 
[22], one can easily obtain the DT invariants for the Calabi-Yau manifolds we use in this 
paper, from [23] and [24]. We will state the DT invariants of interest for the models we 
study, where they are directly applied, in chapters 3 and 4. 

As was done in [1], we will use the index for D4-D2-D0 BPS states of total charge F 
from [3] to enumerate states: 



with the sum running over all possible flrst splits T — t- Pi +r2, belonging to a full split flow 
tree, and (ri,r2) is the symplectic intersection of the two charges, as deflned in appendix 
A. The microscopic logic behind this formula is that all degrees of freedom in a D6/D6 
can be factorized as the degrees of freedom on the gauge theories of the D6 and D6 plus 
the degrees of freedom of the tachyon fleld, which are counted by the intersection product. 
This formula is not accurate though, for all cases, as the moduli space of a bound state is 
in general a fibration: this is discussed in detail, in chapter 4. 

In general, as discussed in appendix C, a specific charge can give rise to several split 
flow trees. A split flow tree will contribute a term to the index of the D4 system as follows: 




(2.1) 



r^ri+ra 



AniTDi) = (-l)l<rD6,ri56>|-i 



I {Tdg, Pije)! iVDT(/3i, ni) A^dt(/32, ^2) 
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where 

n = \x{Cp)+N (2.3) 

with n the total and the added DO-charge. 

The research presented in this paper is closely connected to an 'OSV-like statement' for 
D-particles. The topological string partition function is used to enumerate D6 (and lower 
dimensional brane charge) systems, or, in the mirror picture, to count D3-brane systems. 
Schematically, the relationships, of which we would like to gain a better understanding, 
can be expressed as 

^i'particle ~ I ^top ^ ~ | ^^DT ^ ~ | 2^1)6 | ^ = \Z D?,?' , (2.4) 

where: 

• 2^_Dp^rticic ~ l^topP is an OSV-like statement for a D-particle. 

• I ^top P ~ I -^DT I ^ expressed that one uses the asymptotic expansion of the topological 
string partition function (which is not known exactly) known as the Donaldson- 
Thomas (DT) partition function. 

• |2^DTp ~ |^D6p expressed that one uses DT invariants to enumerate D6-D4-D2-D0 
BPS states^ 

• l-Z^jgp = |^D3p expressed that mirror symmetry comes into play, and that we ex- 
amine the mirror D3-brane systems, as these allow computation of exact central 
charges. 

The partition function of a D-particle ^^Dp^^rtidc identical with the modified elliptic 
genus, associated to the M-theory picture of the worldvolume description of these ensembles 
of BPS D-brane states. Elliptic genera are roughly sketched in appendix B. The rough 
scheme discussed above is followed up in this paper, resulting in various elliptic genera for 
specific CY 3-fold study models, obtained from a new perspective. 

2.2 Polar states: the chromosomes of D particles / black holes 

For partition functions of D-particles (whether they have enough mass to backreact a black 
hole or not) modeled using a mixed ensemble of branes^, the degeneracy of a finite number 
of BPS microstates, the polar states^ determines all other degeneracies through modularity. 
This is illustrated in figure 1. 

Although there are of course many shortcomings in the analogy we are about to pro- 
pose, there is a similarity between how the degeneracy of polar states dictates 'the whole 
rest of possible microstates of a D-particle', and how knowledge of the genome on chromo- 
somes in microbiology dictates — we allow ourselves to simplify things considerably — 'all 

^The closer relationship between D6-D2-D0 states and DT invariants has been clarified in [3]. 

^In the original OSV paper, [11], the authors argue, that the choice of a mixed ensemble is natural from 
the viewpoint of topological string theory (TST), as electric vs. magnetic for topological strings is like 
a choice of position vs. momentum. TST computes topological invariants enumerating constituent BPS 
states, thus there is also some motivation for the connection of TST with BPS bound states. 



-7- 



random BPS microstate 



polar state 



Figure 1: A D-particle's entropy is explained by the many possible ways to model BPS microstates 
by wrapping branes around cycles in compact dimensions. Among the many BPS microstates of 
a D-particle modeled as a mixed ensemble, a finite number of polar states (pictured on the right) 
determine the entire partition function. Polar states come as bound states: the two branes (red 
and blue) are glued together by tachyonic string modes (green). 

possible molecular states of an organism'. We thus like to think of polar states as a sort 
of 'chromosomes for D-particles (and black holes)'. Additionally, just like chromosomes 
come in pairs for (nearly all) mammals ('diploid' organisms), polar states are realized as 
bound states, thus they come in pairs of 'parent chromosomes'. This is illustrated in the 
following figure: 



Figure 2: Polar states as chromosomes: The finite number of polar states determine the entire 
partition function of a D-particle/black hole, and they appear as pairs (held together by tachyonic 
string modes). Also chromosomes for diploid organisms come in pairs, and they encode the entire 
information of the 'microstructurc' of the organism. 

We will come back to this metaphor, when we discuss our results in section 6. Of 
course, these thoughts are nothing else than some intuitive interpretations of the modular 
properties and the information encoded by the poles of a weak Jacobi form. For now, 
this may serve as a motive to focus on a D-particIe's polar states. We will however also 
investigate some of the 'almost polar' states^, which will provide the most interesting checks 
on our methods. 



*By almost polar, we mean non-polar states, very close to the polarity-bound, thus carrying charges 
which are very close to charges of polar states. 
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3. Elliptic genera from split flow trees and DT invariants 



This section is a short account of our results on all polar states for two Calabi-Yau (CY) 
varieties, which allows us to predict the corresponding elliptic genera. We start with a CY 
given as a sextic hypersurface. 



3.1 Polar states on the sextic hypersurface in VFIPnii2 

In VFP|^^]^2; the adjunction formula shows that one can obtain a CY hypersurface by 
choosing a degree six polynomial. The total Chern class reads c{X) = ^^^^^ (i+2Jj") _ 
1 + 14ij2 _ Q8H^, and using = fwipf^^^^ = ^' g^*^ x{X) 

J-^ iJ^ = 3, the weak Jacobi form is in this case three-dimensional: 



i+m 
-- -204. 



As 



Z{q, q,z) =J2 ^^'^^^ ' 



(3.1) 



k=0 



which means that we only have to determine Zq and Zi . By choosing a basis for Lj^ one can 
easily see that there are two gluing vectors, but by the symmetry 7 = —7 one knows that 
one has just one 'fundamental' gluing vector 71. According to our previous intuition, this 
means that in order to determine the complete elliptic genus, one will have to enumerate 
states in the classes [0,%] and [71,(70]- For convenience, we list the DT invariants for the 
sextic of interest: 



Donaldson Thomas invariants: sextic 




n = 


n = 1 


n = 2 


n = 3 


p = o 


1 


204 


20'298 


1'311'584 


/3 = 1 





7884 


1'592'568 


156'836'412 


/3 = 2 


7884 


7'636'788 


1'408'851'522 


136'479'465'324 


P = 3 


169'502'712 


443'151'185'260 


5'487'789'706'776 


440'554'251'409'968 



1. Aq=0,Aqo = 0, 



fO 



45] 



The pure D4-brane carries half a unit of flux to ensure anomaly cancellation, and 



has total charge T = H + ^ + { 
we have introduced the notation 



x{P) 



24 



I UJ 



H+lH^ + lH^ 



(0,1,|,|), where 



(3.2) 



which we will equally use from now on. As in [1], we will also label charge systems by 
their deviation in D2-brane charge Aq and DO-brane charge Ag'o as measured from 
the most polar state. In the 'charge shift' notation it is denoted as Aq = 0, Ago = 0. 
As explained in [3], various charges are related by flux shifts. Charge equivalence 
classes contain the same entropy, and they can be labeled by the (flux) gluing vector 
(see appendix B), as well as the reduced DO-brane charge. A charge equivalence class 
can thus be labeled by [7, go]) and the most polar state for the sextic lies in the class 
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[0, ||]. One finds a split flow tree with centers 



^ 13 9. 

^ ' ' 4 ' 4' 



r2 = (-1,0,-^,0). 



It looks and is enumerated as follows: 



D6h D6 



-4. 



The BPS index reads 

n = (-l)l<ri-r2>|-i|^ri,r2)|iVDT(0,0) • iVDT(0,0) = • 4 • l • l = -4. (3.3) 

Note that the fact that the intersection number between Fi and r2 nicely corresponds 
with the fact that the moduli space of the hyperplane if C X is a CP'^, hence 
;\;(CP^) = 4, because the coordinate with weight 2 can of course not be used to define 
a hyperplane. 

2. Aq=0,Aqo = -l, [0,i]: 

Adding one DO, one gets the total charge (0, 1, 1, |), with reduced DO-brane charge 
go = Ij- The flow tree is analogous to what was found for the quintic [l](the side of 
the DO after the first split can be chosen, according to where one is with respect to 
the appropriate threshold wall) . The charges of the centers after the first split read 



= (-1,0,-^,-1) 



and the flow tree looks like 




DO 



612, 



D6h D6 

where of course 

n = (-i)i<ri.r2>i-i|^ri,r2)|iVDT(o, o) • iVDT(o, i) = (-1)' • 3 • 1 • 204 = 612. (3.4) 

3. Aq=l,Aqo = -l, [71,^]: 

One can now consider a flux, which will involve what we called the relevant gluing 
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vector 7i. According to our intuition, which receives further support at this point, 
this means turning on an extra flux dual to a degree one rational curve. This leads to 
the total charge (0, 1, |, |), and to the reduced DO-brane charge % = thus, there 
is only one polar state in this 71-class. One finds the split flow tree with a pure D6 
one one side, and a D6 with a D2 on a degree one rational curve, as expected. The 
charges read 

r: = (l,l,if.^), 



with the split flow tree 



= (-1,0,-^,-1), 



T = = -15'768. 

/v 

The BPS index is calculated according to 

O = (-l)l<^i'^2>|-i|^ri, r2)|A^DT(0, 0) • Nbt{1, 1) = (-1)^ • 2 • 1 • 7'884 = -15'768. 

(3.5) 

Using a basis for modular forms of the right weight, one can use these numbers to determine 
the modular form to be given by 

Zo{q)= g-S(-4 + 612g-40'392g2 + 146'464'860g^...) (3.6) 

Zi{q) = Z2{q)= g-i(-15'768g + 7'62l'020g2 + ...)), (3.7) 

whose uniqueness follows from [3, 13, 14]. This agrees with the findings of [2] (up to an 
overall sign), which is of course not a surprise, given that the small number of polar states 
supporting split flow tree realizations apparently do not involve subtleties. 

3.2 Polar states on the octic hypersurface in W^lPfiii4 

For iyP|^;^]^4, the adjunction formula shows that one can obtain a CY hypersurface by 
choosing a degree eight polynomial. The total Chern class reads c{X) = ^"^"'"'^'^g^^''^^ = 
1 + 22i?2 _ i48i/3^ and using i?^ = /^yp4^^^^ 8H^ = 2, one gets xiX) = -296. As 
J-^ = 2, the weak Jacobi form is in this case two-dimensional: 

1 

Z{q,q,z)=Y,Zkiq)@k{q,z), (3.8) 

k=0 

which means that we have to determine Zq and Zi. By choosing a basis for Lj^ one can 
easily see that there is only one gluing vector, 71. According to our previous intutition, this 
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means that in order to determine the complete ehiptic genus, one will have to enumerate 
states in the classes [0,%] and [71,^0]- For convenience, we also list the DT invariants of 
interest, for the octic: 



Donaldson— Thomas invariants: octic 




n = 


n = 1 


n = 2 


n = 3 


^ = 


1 


296 


43'068 


4'104'336 


13 = 1 





29'504 


8'674'176 


1'253'300'416 


13 = 2 


564'332 


204'456'696 


45'540'821'914 


6'127'608'486'208 


/3 = 3 


8'775'447'296 


6'313'618'655'104 


1'225'699'503'521'536 


141'978'726'005'461'504 



We will be quite brief on the summary of results on polar states, given the analogy to 
the previously discussed cases. 



1. Aq = 0,Aqo = 0, 



10 



23] 
' 12 J 



The most polar state is the D4-brane carrying flux y for anomaly cancellation, 
with total charge (0, 1, 1, ^). The reduced DO-branc charge can be calulcated to be 



qo - 12 
centers 



and the state is of course in the class [0, ||]. One finds a split flow tree with 



with split flow tree 



r2 = (-1,0,-^,0), 



2 



D6h 



D6 



where the BPS index is calculated as 

= (_l)l(ri,r2>|-i|^r^^r2^|^j^^(0,0) •iVDT(0,0) = (-l)=^-4 - 1-1 



(3.9) 



Note again the fact that the intersection number between Fi and r2 equals —4 cor- 
responds with the fact that one cannot use one of the five coordinates to define a 
hyperplane, in the D4-picture. 



2. Aq = 0,Aqo = -1, 



[0, 



12J 



Adding one DO, one gets the total charge (0, 1, 1, g), with reduced DO-brane charge 
go = X5 • The flow tree is again analogous to what we found for the sextic (the side of 
the DO after the first split can be chosen, according to where one is with respect to 
the appropriate threshold wall). The charges of the centers after the first split read 

17 13 

11 
~6' 



(1,1, 
(-1,0, 



-1 
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and the flow tree looks like 




where of course 



n = (-i)i<'^i'r2)hi|(ri,r2)|iVDT(o,o) • a^dt(o,i) = (-1)^ • 3 • 1 • 296 = 888. (3.10) 

3. Aq= l,Aqo = -1, [71, |]: 

One can now consider a flux, involving the relevant gluing vector 71. As previously 
seen, this means turning on an extra flux dual to a degree one rational curve. This 
leads to the total charge (0,1,2, |), and to the reduced DO-brane charge qq = g: 
thus, there is again only one polar state in this 7i"class. One finds the split flow tree 
with a pure D6 one one side, and a D6 with a D2 on a degree one rational curve, as 
expected. The charges read 

r2 

with the split flow tree 
The BPS index is calculated according to 

n = (-l)l<ri'r2>l-i|(ri,r2)|iVDT(0,0) • iVDT(l, l) = (-1)^ • 2 • 1 • 29'504 = -59'008. 

(3.11) 

Using a basis for modular forms of the right weight, one can use these numbers to determine 
the modular form to be given by 

Zo{t) = g-T5 (-4 + 888g - 86'140^2 + 13l'940'l36g^..) (3.12) 

Zi(r) = 2'2(r) = g-i(-59'008g + 8'615'l68g2 + ...) (3.13) 




This again agrees with the findings of [2] (up to an overall sign). 
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4. Refined enumeration of bound states and elliptic genera revisited 



Using the BPS index for bound states explained in part 2.1, we were able to predict the 
elliptic genera for two study models, and we found exact agreement with the authors of 
[2]. In this section, we will discuss the non-polar states that demand the use of a refined 
prescription for the computation of the BPS index for bound states. In [1], this was 
demonstrated for a non-polar BPS state on the quintic 3-fold: the refined computation 
yielded an exact match with the result predicted from modularity. In our case, the reason 
for the necessity for refined calculations is different, but we again find exact agreement 
with the result predicted by modularity, for all tractable examples. We will show that 
refinements can also alter the enumeration of polar states, in the next section. 

We start with a subsection explaining the main reason why refinement is needed, 
followed by a subsection explaining the more specific techniques used in our computations. 
We then move on to the various explicit examples, and we will also show how our refined 
computations allow us to introduce partitions of the Donaldson-Thomas invariants used 
to enumerate the BPS states belonging to a D6-D2-D0 constituent of a bound state. 

4.1 The main idea: non trivially fibered moduli spaces of bound states 

In the microscopic D-brane picture, a bound state of a D6-brane system carrying a flux 
Fi and bound to N DO's, and a D6 carrying a fiux F2 is held together by tachyonic string 
modes. For the pure D6-D6 bound state, the tachyon field T ([25]) can be understood as a 
map r : — )• F2, where Fi^2 are the line bundles corresponding the fluxes. Alternatively, 
the tachyon can be understood section of 



and the Riemann-Roch theorem allows us to compute an index Xt, counting the number 
of basis elements of this space of sections, r(F| (8) Fi): 



where ch.{F) denotes the Chern class of F, and Td{X) denotes the Todd class for the base 
space (in this case the CY). If one adds a DO-brane, the tachyon field needs to vanish at an 
additional point. The map is then a section of T{F2 (S> Fi x Xp), where Ip denotes the ideal 
sheaf on the point where the DO resides. When N DO's are bound to the D6~brane, the 
tachyon field T needs to vanish also at these points, which means that one has to impose 
N constraints on T: the number of independent sections of the tachyon fields is reduced 
by (as opposed to the case where one studies a D6-D6-system without any extra bound 
DO's). In this way an index is computed, accounting for the degrees of freedom associated 
with the tachyon moduli space A^t- 

The index presented in subsection 2.1 is based on the assumption, that the moduli 
space of a bound state factorizes into the tachyon moduli space, the moduli space of the 
D6-brane, and the moduli space of the D6-brane: 



Ter{F^(g)Fi) 



(4.1) 




(4.2) 



(4.3) 
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To compute an index one naively uses, [3] (compare also with subsection 2.1), 
(_l)l(ri,r2>||^ri,r2)| • 7VDT(/3i,ni) • A^dt(/32, na), 

where (— l)l^'"i''"2)l|(ri,r2)| and Xt agree up to a sign. 

In general, this moduli space is non-trivially fibered [1], and thus the dimensionality 
of the fiber can perform jumps. As a consequence, the index 4.2 is not always accurate: 
it computes a virtual dimension [26], which does not always equal the real dimension of 
the moduli space. Basically the reason is that the tachyon fields do not perceive all the 
constituent states generically. For simplicity, we will restrict the discussion in the following 
to the case that only one of the two constituents is not generically perceived, namely 
the D6-consistuent. We will also focus on the case where DO-branes are not perceived 
generically by the tachyon fields, although we will show for an example, that an analogous 
phenomenon holds for special D2-states (and the curves on which these are wrapped). 
Generalizations of the presented scheme will become clear from our examples. 

We show that the moduli space of a BPS bound state splits up for our examples, and 
can be grouped into two pieces: 

The first part, with the superscripts 'g', stands for the part where the tachyon fields perceive 
the constituent states 'generically', and hence the virtual dimension of A^t is actually the 
real dimension, and the part with the superscripts 's' stands for the part where the tachyon 
fields perceive the constituent states as 'special states' (in the present case, always for the 
D6-system). This usually happens because the tachyon fields do not perceive all DO's that 
are bound to the D6, which results in a constraint loss on the tachyon field T. For example, 
generically three DO's mean that the independent sections have to vanish at three points. 
If the tachyon is blind to one of the three, this number is reduced to two, resulting in a 
jump in the dimension of the fiber of the moduli space of the bound state. Typically, the 
dimension of Mt^ will thus be greater by one, as opposed to A^Tg- 

In [1], it was proposed to define Donaldson-Thomas densities (in analogy to a top 
Chern class), that integrate over moduli space to a Donaldson-Thomas invariant, in order 
to be able to define a product formula for the bound state index on the level of index 
densities. One can also decide to integrate these densities on the various partitions of 
moduli space that have a constant dimension of the fiber of the tachyon moduli space. 
We call these integrated densities Donalds on- Thomas partitions (DT partitions). They 
separately enumerate generic and special D6-D2-D0 states, and we will calculate a series 
of examples of such DT partitions. Using these, one can write down how our refined indices 
come about: 
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(_l)l(ri,r2>||^ri,r2)| • iVDT(/3i,ni) 



Mt^ Md& 




A/'S(^2,n2). 





Before giving our explicit examples, we will now give a more detailed instruction on how 
to separate the 'generic' from the 'special' states, resulting in a calculation of Donaldson- 
Thomas partitions. 

4.2 Algebraic techniques to deal vi^ith special constituent states 

In [1], D4: — 3D0 states had to be treated more carefully, because the tachyon fields of a 
D6-(D6-3D0)-bound states do not perceive three collinear DO's generically. Rather, they 
appear only as two particles. In terms of algebraic geometry, this is simple to express: 
the three constraints on the tachyon field are not independent (only two of them are 
independent). For our examples, this constraint loss occurs for two reasons, which we shall 
refer to as 'special loci', and as 'special tangent directions' (which are important when 
blowups are performed, required for dealing with coincident loci for two DO's). The special 
loci occur because the Calabi-Yau varieties are embedded in weighted projective spaces, 
which means that there is not a complete democracy amongst coordinates: If a DO sits at 
a position with non-zero coordinates of a higher weight only, it will not impose constraints 
on the tachyon, as the higher weight coordinates cannot be included in the definition of 
the tachyon map. Let us elaborate on this in more detail. 

1. Constraint loss because of special loci 

For simplicity, consider a Calabi-Yau embedded in VFIPini^ with n > 1, coordinates 
(xi,...,X5) and a bound state of a D6-brane with flux F = H (one unit of flux) 
and a D6 without flux. In this case, H \s a. line bundle of which the coordinates Xi 
form sections. However, X5 is forbidden as a section, as it has a higher weight. This 
means that there are only four instead of five independent sections for the tachyon 
T G r(i/) (F| is the trivial bundle in this case). The most general tachyon field 
reads T = aixi + ... + 04X4. In general, placing a DO on the D6 means imposing one 
constraint on the tachyon field. A good way to think about this, is by treating the 
DO on the D6 with point particle quantum mechanics, (check [27]). If one puts the 
particle at X5 = 1 (and all other coordinates zero), it will not impose a constraint on 
the tachyon field. This might be of relevance and it might not: one still has to check 
whether this point actually lies on the Calabi-Yau. We will encounter cases where 
situations analogous to this fictitious example arise. 

2. Constraint loss because of special tangent directions 

This was not the whole story, though. An additional complication arises as soon as 
one considers two (or more) DO-particles: orbifold singularities arise when particles 
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meet. This is dealt with by performing a blowup, which imposes a distinction be- 
tween the particles, that inhabit the same spot on the brane (and hence on the CY). 
One can intuitively picture this as considering an infinitesimally short time period 
before they meet, and distinguishing the particles upon all the different tangent di- 
rections (which of course have to be tangent directions to the Calabi-Yau variety 
under consideration), from which the two particles can approach each other. Two 
particles lying on the same spot in the Calabi-Yau would originally only impose one 
constraint on the tachyon field, according to the previous discussion. After perform- 
ing a blowup, an additional constraint arises, from the splitting through the tangent 
direction. This means, that the tachyon moduli space fiber does not jump (at least 
generically) for states, where particles lie at the same locus. The point is: it can 
jump. Namely, one (or several) tangent directions might be built from coordinates 
that do not impose constraints. Thus, also after performing blowups, one still has to 
distinguish between generic and special states of particles, that lie at coincident loci. 

Let us summarize the two situations, when special states occur (in this case for Do's), once 
more: 

• For non-coincident DO's, one needs to check whether the particles lie at 'special loci', 
where they do not impose a constraint on the tachyon fields. One could refer to these 
cases as the 'special non-blowup loci'. 

• For a bound state, there are cases, when various DO's are coincident: in this case 
one needs to perform blowups. These blown-up states have to be separated into 
generic and special states, according to whether the tangent direction (arising from 
the blowup) imposes a constraint on the tachyon fields or not. One might refer to 
these latter cases as the 'special blowup loci'. 

As will become clear from an example later on, analogous implications arise for bound 
D2-branes, depending on whether the curves (on which those branes are placed) impose 
'enough' constraints on the tachyon fields (to make the bound state 'generic'). 

4.3 Special points and DT partitions 7V^^^*^(0,2) 

In the following, a few interesting non-polar states on the sextic and the octic CY's for 
which the elliptic genera were predicted in the previous section, are examined. The refined 
caculations match the predictions from modularity. 

The state Aq = 0, Aqo = —2 on the sextic 

For this charge system, one finds a split flow tree with centers 

ri = (M,^,^), 
r2 = (-1,0,-^,-2), 
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of the form 




= -2 • 20'298 = -40'596, 



and obviously, the index obtained naively r^naive differs from the exact index, which 
is given by riexact = — 40'392 (see equation 3.6). This can be put right using the refined 
index, as the non-trivially fibered moduli space of this bound state dictates. 

One can either argue from the D6-picture, or from a purely algebraic geometrical 
D4-perspective. Essentially the two arguments are identical. In the first picture, one 
chooses two point-like sheaves to lie on the vanishing locus of the polynomial describing 
the tachyon; in the second case, one places two DO-branes on the D4. 

As dicussed before, for the most polar state, the most general tachyon map is of the 
form 

T = aixi + ...04x4, (4.5) 

choosing coordinates (xi, X4, X5) on VFIP11112 transforming with weights (1,1,1,1,2). 
Setting two DO-branes on the D6 means that the map has to 'pass' through two points, so 
in general, the moduli space of the tachyon will be reduced to from CP'^ to CP^, yielding 
x(CP^) = 2. This is where the intersection number 2 comes from. It will now be shown 
that the dimension of this fiber jumps (as was observed for a state on the quintic 3-fold, 
though for a different reason, in [1]). One might say, that two DO~particles behave as one 
on a specific 'special' locus. The scheme is to analyze the remaining moduli describing this 
tachyon, demanding that two points as well as the tangent directions at these points lie on 
the CY as well as on the zero locus of the tachyon. 

Finding special tangent directions arising from blowups 

We discussed in the previous section that the fiber of the tachyon moduli space jumps 
in two cases. Firstly, if a particle is placed at a 'special locus', directly. Secondly, if two 
particles are coincident and the tangent direction arising from the blowup procedure is a 
'special tangent direction'. Before putting all the bits into place for the state under current 
investigation, we will look for the special tangent directions for the blowups which we have 
to perform. Before continuing, let us stress that this is a mechanism special to Calabi-Yau 
varieties embedded in weighted projective spaces, but will probably play an even more 
prominent role when working with more general Calabi-Yau varieties, embedded in toric 
ambient spaces. 

Generically, demanding that two points lie in the zero of (4.5), imposes two independent 
constraints, reducing the number of moduli by two. Denote the coordinates of the two 
points by x^^ and x^^. This is not that clear if the two points lie on top of each other 
x^^ = x^'^ . However, this is generically resolved by the blowup procedure, where the 
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tangent direction of the blowup gives an independent constraint, again leading to two 
constraints. However, as one is working with a weighted projective space, the coordinate 
X5 does not appear in (4.5). Different xs-coordinates are not 'seen' by the tachyon. The 
points do however have to lie on the sextic, for which one can choose a representative given 
by a transverse polynomial, [28] , 

Psextic = +p('^^(xi,X2,X3,2;4) = 0. (4.6) 

where p*-^^ is a degree 6 homogeneous polynomial in the indicated coordinates. This means 
that the case x^^ = e 3 for k = 0,1,2 with all other coordinates equal (xf^ = x^^ for 
i = 1, 2, 3, 4), is possibly of interest. 

If the first four coordinates of the two points are identical, this amounts to imposing 
only one instead of two constraints on the tachyon field. After the usual blowup procedure 
(if the fifth coordinate is also equal), the particles are however distinguished by a tangent 
direction, which can be interpreted as the direction from which they 'approach' each other. 

Consider thus a tangent vector X^di, and demand that it is a direction tangent to the 
sextic CY; 

VxPsextic = 0, (4.7) 

But it also acts on (4.5) 

Vxr = X^ai + ... + X^a4. (4.8) 

Setting S/xT = 0, is like a second constraint. The particles are distinguished by the 
tachyon field after a blowup. This is also the reason why a coincident locus does not in 
general lead to a jump in the tachyon fiber above that locus. Namely, (4.8) does not lead 
to an extra constraint, iff 

/ XI X2 X3 X4 \ , . 

rank(^^,^,^3^,j<2, (4.9) 

which can happen only if either = ... = X^ = or if X* = Axj, with A G C*. We will 
show where this can happen in the present case. In the following distinction of cases, as 
well as in similar computations later on in this paper, we will always refer to coordinates 
of the coincident particles. X5 7^ for example means, that the two (a priori) coincident 
particles lie on a locus with X5 7^ 0. In the present case, one can distinguish between 

• X5 7^ 0: This means one can choose affine coordinates with X5 = 1. Thus, in these 
coordinates, one knows that X^ = for the tangent vector (and hence the case 
X^ = ... = X'^ = is ruled out and only the case X* = Axj remains). The tangent 
vectors should remain in the CY (4.7), which leads to 6Xp'^^^ = 0. Since A 7^ 0, 
this means p^^^ = 0, so X5 = upon plugging this into (4.6), which contradicts our 
assumption. If X5 7^ 0, the two particles are distinguished after the blowup. 

• X5 = 0: In this case, fix (w.l.o.g.) xi = 1. This means X^ = for the tangent 
vector. To have only one instead of two 'constraint' equations for the tachyon, one 
thus needs X^ = ... = X^ = 0. This occurs when the tangent vector equals X^d^, 



-19- 



which is possible at the locus X5 = 0. Note that there is one tangent direction and 
one locus (X5 = 0) for which this happens. This is where the fiber of the tachyon 
moduli space jumps and needs to be taken into account. 

Calculating the refined index 

Starting from the naive form for the index of this bound state, 

x(two particles) • x(tachyon) = xi^)"^ " x(CP^) + corrections, (4-10) 

one can nicely see how the refinement comes in. The refined index receives the following 
contributions (where x{^) = —204 denotes the Euler character of the sextic). 

• ^{x{X)'^ — 3x{X) + 2xo)-x{'^^^)- this is the generic locus, where the two particles are 
separated and the locus where xf^ = x^^ for i = 1, 2, 3, 4 has been substracted. Note 
that one has to be careful not to substract the locus where the first four coordinates 
are identical and X5 = more than once. This has been taken into account with 
the +2x0 term, xo is the Euler character of the locus X5 = in the CY and can be 
calculated from the adjunction formula to give xo = 108. The factor ^ accounts for 
the fact that interchanging the two particles gives the same configuration. 

• (x ~ Xo) • xlCIP^) • xlCIP^): this accounts for the locus where the two DOs coincide, 
without the locus = 0. Note that the x(CP^) results from the blowup of a 
codimension 3 locus. 

• 2 • ^(x — Xo) ■ x(dP^): this takes into account the loci where x^^ = e^s^x^^ for A; = 1 
and k = 2 (hence the overall factor of two, as these loci both contribute equally). 
Note that the fiber of the tachyon has jumped, these two particles are seen as one. 

• Xo • (x(CIP^) - 1) • X(CP^): here, the locus xf^ = xf for i = 1,2,3,4 and xf = is 
dealt with. In principle, one just has to do a blowup of a codimension 3 locus (hence 
a factor of x(CP^)). After the blowup, the tachyon 'sees' two DO's. However, one 
needs to substract the one tangent direction we found in the above analysis, because 
one loses one of the two 'constraints' on the tachyon. This one tangent direction is 
taken into account on the next line. 

• Xo • 1 • x(CP^): for this blowup direction (for which the 1 stands for the index), the 
tachyon again sees only one particle. 

Note that had one ignored the subtlety with the locus X5, one would have retrieved the 
calculation where one finds a product of the tachyon index 2, and the Donaldson-Thomas 
invariant iVDT(0, 2). Collecting all the pieces linked to the value 2 or 3 for the tachyon 
index (up to a sign), one can state the correct index in the form 

^^exact = -2 • (20'298 + 204) - 3 • (-204) = -40'392. (4.11) 
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This allows us to state the Donaldson-Thomas partitions Afjyr i^^'^) sextic: 

A/'i?J(0,2) = 20'502, (4.12) 
A/'^"^(0,2) = -204. (4.13) 

J\f^^{0, 2) counts the generic D6-2D0 BPS states, for which the tachyon perceives two DO's, 
and Mj^r^[0, 2) counts the special D6-2D0 BPS states, for which the tachyon perceives only 
one DO-brane. Note that there is a sign difference between these indices. 

ArDT(0,2) =<?)(0,2) +<^^(0,2). (4.14) 

Let us consider a second example. 

The state Aq = 0, Aqo = —2 on the octic 

For this charge system, one finds a split flow tree with centers 



of the form 




= -2 • 43'068 = -86'136. 



For the octic, one can choose 

Poetic = x| +p(^)(xi,X2,X3,X4) = (4.15) 

as a transverse polynomial. Again, the cases where the particles have four equal coordi- 
nates, and where x^^ = e 2 , are possibly of interest. The analysis of the locus where 
the tachyon field perceives the BPS state with two DO's differently and where the fiber 
changes dimension is analogous to the sextic case, and the discussion will therefore be 
brief. We again start by searching for the 'special tangent directions' by considering (4.9): 

• X5 7^ 0: This means one can choose affine coordinates with 2:5 = 1. Again, one knows 
that = for the tangent vector (and hence the case = ... = X'^ = is ruled 
out and only the case X^ = Xxi remains). The tangent vectors should remain in the 
CY, which leads to 8Xp^^^ = 0. p^^^ = would imply 3:5 = 0, contradicting the first 
assumption. 

• X5 = 0: In this case, fix (w.l.o.g.) xi = 1. This means X^ = for the tangent 
vector. To have only one instead of two 'constraint' equations for the tachyon, one 



- 21 - 



thus needs = ... = = 0. This occurs when the tangent vector equals X^d^, 
which is possible at the locus X5 = 0. Note that there is again one tangent direction 
for which this occurs. This is where the fiber of the tachyon jumps and needs to be 
taken into account. 

Calculating the exact index 

The index receives similar contributions as in the case of the sextic, but the calculation 
is slightly simpler. 

• ^{x{Xy — 2x{X) + xo) ■ x(CP^): this again is the generic locus, but as in this case 
x^^ = e~2~x^^ for A; = 0, 1 one substracts two instead of three loci with index x{^)- 
Instead of substracting the locus 2:5 = three times, one does this twice, and needs 
to compensate once. In this case, xo = 304. 

• (x — Xo) ■ x(dP^) • x{C^^)- this accounts for the locus where the two DOs coincide, 
without the locus = 0. Note that the x(dP^) results from the blowup of a 
codimension 3 locus. 

• ^(x ~ Xo) ■ x(dP^): this takes into account the locus where x^^ = —x^^ and other 
coordinates equal. This is again a locus where the tachyon index has jumped. 

• Xo • (x(dP^) — 1) • x(dP^): here, the locus X5 = and other coordinates equal is 
dealt with. In principle one just has to do a blowup (hence a factor of x(CP^)). 
After the blowup, the tachyon 'sees' two DO's. Again, one needs to substract the one 
tangent direction we found above, because one loses one of the two 'constraints' on 
the tachyon. This tangent direction is taken into account on the next line. 

• xo-l-x(CP^): for this one blowup direction, the tachyon again sees only one particle. 

Collecting all the pieces linked to the value 2 or 3 for the tachyon index (up to a sign), one 
can state the correct index in the form (to be compared with equation 3.12) 

f^cxact = -2 • (43'068 - 4) - 3 • (+4) = -86'140. (4.16) 

This means that the Donaldson-Thomas partitions J\fj^r^\o,2) for the octic read: 

A/'^?)(0,2) = 43'064, (4.17) 
a4^^(0,2)=4. (4.18) 

Note that the sum of the partitions yields iVDT(0, 2). 

4.4 Special curves and DT partitions J\f^rp\l,2) 

Up until now, refinements were presented, which were necesary due to the fact that the 
tachyon fields did not perceive DO's generically. An example will be given now, where an 
analogous refinement is necessary because of the fact that a D2/D0-state (with the D2 
wrapped on a curve) is not perceived generically. A different way of expressing this is to 
state that there are special D6-D2/D0 bound states. 
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The state Aq = 1, Aqo = —2 on the octic 



For this charge system, one finds a spht fiow tree with centers 

r2 = (-1,0,-^,-2), 

of the form 




8'674'176. 

D6h m- D2{Ci3=^),m 

with an index naively calculated as 

n = (-l)l<^i'^2>|-i|^p^^r2)|iVDT(0,0) • iVDT(l,2) = (-1)° • l • l • 8'674'176 = 8'674'l76. 

(4.19) 

This naive index agains needs refinement. In this case, for simplicity, choose the 
Fermat-polynomial for the octic: 

Poetic = xl + xl + xl + xl + xl = Q (4.20) 

The following discussion is nonetheless also valid for more general transverse polynomials. 
A degree one rational curve on the octic can be represented as a degree one map from CP^ 
to the Calabi-Yau. Consider for example the map 

(s, t) (s, eTs, t, e^t, 0). (4.21) 

This generically imposes two constraints on the tachyon field, reducing its moduli 
space to CP^. Adding an extra D-particle will then reduce this moduli space to CP*^, 
unless something special happens: 

• The particle (DO) does not sit on the curve, but nevertheless produces no extra 
constraint. It is easy to verify that this cannot possibly happen for this example. 

• The particle lies on the curve, which means that a blowup needs to be performed 
in the directions normal to the curve. Again, one might encounter special tangent 
directions, which do not impose an extra constraint on the tachyon. Following a 
similar procedure as in the previous examples, one can indeed verify that this is the 
case for the direction X^d^. As 2:5 = lies on the curve (4.21), this direction is 
automatically also tangent to the octic. 

Calculating the exact index 

The various contributions to the exact index according to the refined prescription read: 
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• N£)t{^, ~ Xc)x(dP ), where xc = 2 is the Euler characteristic of the curve. 
This term deals with the case, when the DO is placed at a locus different from the 
curve, thereby reducing the tachyon moduli space to CP'^. 

• 7Vdt(1, 1)xc[x(CP^) - l]x(CP°), dealing with the case that the DO is located on the 
curve, but the blowup tangent direction leads to an extra constraint on the tachyon. 

• N£)t{1, 1)xc ■ 1 ■ xlCP"*^), which deals with the case, when the DO lies on the curve, 
and a blowup is performed leading to a special tangent direction. This is an example 
of what was referred to as a special D2/D0 bound state, or alternatively just as a 
special curve. In this case, the tachyon field moduli space remains a CP^. 

In total, this leads to the index 

^^exact = 1 • (8'674'176 + 59'008) + 2 • (-59'008) = 8'615'168. (4.22) 

Spectacularly, by comparing this number to the prediction from modularity (3.12), one also 
finds exact agreement for this case! One can thus state the Donaldson-Thomas partitions 
Mbt{^, 2) for the octic: 

A/']^f^(l,2) = 8'733'184, (4.23) 
A/^,'^(l,2) = -59'008. (4.24) 

Again, note that the sum of the partitions yields A'dt(I) 2). 

To summarize, we conclude that all our results obtained in this section provide exact 
agreement with the predictions from modularity. This clearly is strong evidence that our 
procedure is correct. These results also provide a non-trivial and successful test for the 
split attractor flow conjecture. We will now discuss states on another CY manifold, for 
which also a polar state requires calculation of a refined bound state index. 



5. Refined predictions for elliptic genera 

In this section, we show how our refined index computation alters the prediction for an 
elliptic genus of a CY, realized as a degree ten hypersurface in VFPiii25- We choose coordi- 
nates (xi, X4, X5) with the weights (1, 1, 1, 2, 5). We will refer to this CY as the decantic, 
in the following. The total Chern class of this space reads c{X) = (^+-^) = 
1 + 34if2 _ 288i?3, and using = J^^^ ^ WH^ = 1, one obtains x(X) = -288. We 

again wrap a D4-brane on the hyperplane class divisor P = H. For the lattice of fluxes, one 
finds in this case that the pullback Lx = i*p{H'^ {X ,1,) and it's orthongonal complement, 
Lx ® Lj^, are already unimodular, thus no gluing vectors exist and the elliptic genus is a 
'one-dimensional vector', 

Z{q,q,z) = Zo{q)eo{q,z) . (5.1) 
Again, the list of DT invariants of interest is stated below. 
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Donaldson— Thomas invariants: decantic 




n = 


n = 1 


n = 2 


n = 3 


p = o 


1 


288 


40752 


3'774'912 


13 = 1 


1150 


435'827 


89'103'872 


11'141'118'264 


P = 2 


-64'916'198 


40'225'290'446 


9'325'643'249'563 


1'119'938'319'168'004 



1. Aq=0,Aqo = 0, [0,|f]: 

As usual, the most polar state is the D4-brane carrying flux ^ for anomaly cancella- 
tion, with total charge (0, 1, g, The reduced DO-brane charge can be calulcated 
to be go = H, thus the state lies in the class [0, ||]. One finds a split flow tree with 
centers 



3, 



(-1)2 • 3 • 1 • 1 = 3. (5.2) 

Again, the number 3 can nicely be understood along the previous lines. In W/'IPiii25 
one can only use three coordinates to define a hyperplane, and x(Cp2) = 3, so 
the correspondence between divisor moduli in the D4— picture and tachyonic degrees 
between the D6 and the D6 again works out. 

2. Aq = 0,Aqo^-l, [0, ||]: 

Adding one DO, one obtains the total charge (0, 1,|,^), with reduced DO branc 
charge Qo = The flow tree is again analogous to the previous findings, the side 
of the DO after the first split being governed by the appropriate threshold wall. The 
charges of the centers after the first split read 



with split fiow tree 



, 23 19, 
= ^^'^'12' 12^' 

r2 = (-1,0,-^,0), 



DAi 



A 

D6h D6 



where the BPS index is calculated according to 

n = (-i)i<ri.r2>i-i|^ri,r2)|iVDT(o,o) •iVDT(o,o) 



The flow tree looks like 




J, D 

k 



-576 



from which one is tempted to conclude 



Q = (-l)l<ri.r2>|-i|^ri, r2)|A^DT(0, 0) • A^dt(0, 1) = (-1)1 • 2 • 1 • 288 = -576. (5.3) 



Using these two polar degeneracies, the elliptic genus would be determined to be 



Zo{q) = g-i (3 - 576g + 27l'704g2 + 206'40l'533g^ + 2l'593'767'647g'^...), (5.4) 



which agrees with the findings in [2]. We will now argue that this is not quite correct. We 
will predict the elliptic genus to be 



Zo{q) = g- 24 (3 _ 575g + 27l'95552 + 206'406'410g^ + 2l'593'817'025g'^...). (5.5) 



The reason for this lies in the fact that the index (5.3) is not correct, because the state 
Aq = 0, Ag'o = 1 has a non~trivially fibered moduli space. This will allow us to calculate 
partitions of A^dt(0, 1), for the decantic. 



Note that the tachyon map for the most polar state is a section of the bundle H, and 
is of the general form 



as the coordinates X4 and X5 are 'forbidden' (weight too high). This yields a moduli space 
with Euler character x(CP^) = 3, accounting for the degeneracy of the most polar state. 
If one adds a DO-brane to the system (this is also discussed in [1] for examples on the 
quintic 3-fold), the tachyon map also has to vanish on an additional point. Inserting this 
in (5.6) generically eliminates one of the moduli, reducing to a tachyon moduli space with 
Euler character x(CPi) = 2. The problem is that one can place the DO-brane at the locus 
xi = X2 = X3 = 0, which is indeed a point lying on the decantic X. Placing the DO on this 
point means that this 'particle' will not imply a constraint on (5.6). The Euler character 
of this locus {xi = X2 = X3) xo can easily be determined: Xo = 1 (this is trivial, as the 
locus is just a point). 

Thus, the correct index for the Aq = 0, Ago = — 1 system reads as follows: 



DT partitions A/^rjf^(0, 1): generic and special D6 — DO states 



T = aixi + ...03x3, 



(5.6) 



Q = -2-{xo- X{X)) - 3 • (-xo) = -575. 



(5.7) 




(5.8) 
(5.9) 
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Again, note that N^TiO, 1) = A/Qrp(0, 1) +A/']5'p(0, 1). The superscript g stands for generic, 
and A/'^rp(0, 1) counts the number of D6-D0 states, which are perceived by the tachyon 
generically. The superscript s stands for special. Accordingly, NQrp(0, 1) counts the num- 
ber of D6-D0 states, where the DO sits at a special locus, where the tachyon matrix does 
not perceive the particle. 

Of course, this also means, that the index for the split flow tree must be stated correctly: 



This leads to the elliptic genus 

Zo{q) = g-it(3 - 575q + 27l'955g^ + 206'406'410g^ + 2l'593'817'025g^..). (5.10) 

It is interesting to note that the authors of [2], after having predicted a slightly deviating 
elliptic genus (as explained above), find 271'952 as a prediction for the number of BPS 
states of the system, which we denote as Aq = 0, Ag'o = —2. This is only off by 3 of the 
modular result we predict, as opposed to the 248 from the result predicted by the 'naive 
elliptic genus' (5.4). This might be seen as an indication that the new prediction is indeed 
correct. Unfortunately, our technique does not allow to check this prediction, as a single 
flow exists for this state and it is thus not possible for us to confirm our elliptic genus with 
absolute certainty, although we do believe that we have collected strong evidence for our 
computational scheme. It is nevertheless interesting to predict Donaldson-Thomas parti- 
tions AA(3.^)(0,2) for the decantic. 

The state Aq = 0, Aqo = —2 on the decantic 

The total charge for this system reads T = (0,1,^,—^), which implies go = ~if- this 
is thus a non-polar state. One finds a split fiow tree with the centers 




-575. 



D6h 



D6 



and a flow tree of the form 

D4h 



2D0 



D6h D6 



2D0 

= 40' 752, 



which would naively yield an index 
n = (-l)l<ri-r2>|-i|^p^^r2)|iVDT(0,0) •iVDT(0,2) = (-1)° • 1 • 1 • 40'752 = 40'752. (5.11) 
For the decantic, one can choose 

Pdecantic = xf + + p^^^^ (xi , X2 , X3) = (5.12) 

as a transverse polynomial. Note that the moduli space for the tachyon was CP^ for the 
most polar state. Generically, this is reduced to CP'^ when placing two DOs, but there are 
a lot of subtleties involved. Namely, the cases when the DOs have three equal coordinates 
xi,X2,x^ and x^^ = e^^x^'^,x^^ = e^^x^^ (with j = 0,1,2,3,4 and A; = 0,1) are of 
special interest. Additionally, the locus xi = X2 = X3 = is special. Again, constraint 
loss will happen on some loci directly, but will also result from blowups, when placing the 
two DOs on the same locus. Thus, we will start by analyzing which tangent directions are 
special, and thus require special treatment when performing a blowup. 
The condition for a contraint loss to occur after a blowup reads 

which can happen either = ... = X"^ = or if X* = Axj. 

• X4 7^ or X5 7^ 0: Assuming that not all coordinates xi,X2,X3 vanish at the same 
time (this case will be dealt with separately), an analysis shows that for the cases 
that either (or both) of the coordinates X4,X5 do not vanish, one finds one tangent 
direction for which a constraint loss occurs. These subcases shall be discussed briefly: 

- X4 / 0,X5 7^ 0: In this case set X4 = 1, thus = 0. One can easily check 
that X^ = X'^ = X^ = is not possible as it would imply X^ = 0. Thus, 
set X* = Xxi {i = 1,2,3), which leads to VxPdecantic = 10Ap(i°) + 2X^x5 = 0. 
Combining this with (5.12) yields X5 + 1 = ^1^X^x5. This completely fixes the 
tangent vector, thus there is one direction for which constraint loss occurs. 

- X4 7^ 0,X5 = 0: Set X4 = 1, thus X^ = 0. Choosing X^ = X^ = X^ = yields 
VxPdocantic = 2X^X5 = 0, thus One tangent direction. 

- X4 = 0, X5 7^ 0: Set xg = 1, thus X^ = 0. Choosing X^ = X^ = X^ = yields 
VjfPdecantic = 5X^x1 = 0, thus One tangent direction. 
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• X4 = 3^5 = 0: In this case, one can choose xi = 1, and thus X = 0. So, a general 
tangent vector reads ^252 + ... + X^d^. Setting X^ = X"^ = X^ = 0, and plugging 
this into VxPdecantic = yields 5X^x4 + 2X^x5 = 0, which is always satisfied. Thus, 
the tangent directions for which there is a constraint loss, form a CP^. 

• xi = X2 = X3 = 0: In this case, set X5 = 1. In that case (5.12) reads x\ + l = Q. One 
might think that one has found five points on the CY, but taking the equivalence 
relation into account under the group action, one realizes that this is only one point. 
In this case X^ = X^ = 0, and thus X^ = X"^ = X^ = is not possible. Therefore, 
there is no (extra) constraint loss when considering a blowup, when the two DOs 
coincide at this point on the decantic. 

Calculating the exact index 

Using the adjunction formula, one can calculate the Euler character associated to a 
number of loci of interest for the following calculation: 

1. X4 = : X4 = 76. 

2. X5 = : X5 = 295. 

3. X4 = = : X45 = -70. 

4. xi = X2 = = Q : Recall that this is only one point. The Euler character is thus 
Xo = 1- 

A careful calculation reveals the following contributions: 

• \{{x{X) - Xo? - 10(x(^) - Xo - X4 - X5 + X45) - 2(x4 - X45) - 5(X5 - X45) 

— X45) • x(dP*^): this is the generic locus. The case xi = X2 = x^ = Q has been 
substracted from the beginning on, and additionally, also the cases when the coor- 

p livij p p lirik p 

dinates xi,X2,xz of the two DOs are identical and 

j = 0,1,2,3, 4; /c = 0,1 will be treated independently. These cases have been sub- 
stracted, but for each possibility (j, k), the subloci where X4 = 0, X5 = or both have 
been substracted and then treated separately. 

• {x~Xo~ X45) ■ (x(dP^) — 1) • x(CP''): this is the most general case when the two DOs 
coincide. One has to treat various loci separately: the case, when the two particles lie 
on the point xi = X2 = X3 = (this is substracted by the term — Xo); and also the case 
when X4 = = has been removed (the term — X45) a-^d will be treated separately. 
The factor x(CP^) arises from the blowup of a codimension 3 locus. According to 
the analysis presented above, one does however need to substract one direction, for 
which there will be a constraint loss. In this case, two constraints on the tachyon are 
imposed, reducing the moduli space to CP°. 

• {x~ Xo~ X45) ■ 1 ■ x(CP"'^): This is the case analogous to the previous, but associated 
to the blowup direction yielding a constraint loss. 
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• X45 • (x(CP^) - X(CP^)) • x(CP°): When the two DOs coincide and ^4 = ^5 = 0, 
a blowup is performed, but the previous analysis revealed two tangent directions 
associated to a constraint loss. That locus will be dealt with, next. In this case, two 
constraints are imposed on the tachyon, yielding a x(CP°). 

• X45 ■ x(dP^) • x{CP^)- This is the case, when the blowup is associated to a constraint 
loss, with the two DOs coincident on a locus with X4 = X5 = 0. 

• Xo'X(CP^)x(CP^): When the two DOs lie on the point xi = X2 = xs = 0, there is only 
one constraint on the tachyon, arising from the tangent directions after the blowup. 
Recall that there is no constraint from placing a particle at this point: this subtlety 
appeared already when considering the state DA — DO on the decantic, previously. 

• (x — Xo) • Xo • x(dP^): This contribution arises, when one DO is placed on the locus 
xi = X2 = a;3 = (this DO will not impose a constraint on the tachyon), and the 
other one somewhere else. 

• ^ • 9(x — Xo ~ X4 ~ X5 + X45) ■ x(dP"^): These are the cases, when the two DOs have 
identical coordinates xi,X2,X3, but differ in at least one of the other coordinates, 

p 2'Kij p p 2ivik Po X 1 1 • 1 • 1 

x^^ = e 5 x^^,x^^ = e 2 Xr,^. In these cases, there is only one constramt on the 
tachyon. The cases when 0:4 = or X5 = will be treated separately, though. 

• 1 • |(X4 ~ X45) ■ xl^P^): This is as the previous case, but additionally X4 = 0. 

• 4 • ^(x5 — X45) ■ xlCIP^): Again, the conditions as previously, but with X5 = 0. 

Collecting all the pieces linked to the value or 1 for the tachyon index (up to a sign) , one 
can state the correct index in the form 

^^exact = 1 • (40'752 + 3'127) + 2 • (-3'l27) = 37'625. (5.14) 
(1 2) 

The Donaldson-Thomas partitions A/^)^^ (0,2) for the decantic thus read 

M^^^{0,2) =43'879 , (5.15) 
A/J^(0,2) = -3'127 . (5.16) 

Again, note that the sum of the partitions yields A'^dt(0, 2). 

Clearly, the number 37'625 is still very far off from the modular prediction 271'955. 
The missing states cannot however be calculated at the moment. Namely, one finds a single 
flow for this charge system, so there is little hope of obtaining the correct index exclusively 
using the methods utilized in this paper. It is left as a problem for future research to 
enumerate the number of BPS states corresponding to this single flow, and (possibly) find 
and enumerate other split flow trees. 
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6. Discussion 



In this work, we studied examples of BPS degeneracies of D~particle microstates, arising 
from a mixed ensemble of D-branes wrapped on algebraic one-modulus Calabi-Yau va- 
rieties. We have (presumably) found exact results, directly, and have explained and put 
to use a refined computational scheme, which allowed us also to find exact results for the 
other tractable cases. This indeed provides strong evidence for the split attractor flow tree 
conjecture, stating that (single and) split flows are an accurate classification of the BPS 
spectrum of type II string theory. 

We will discuss whether one might be able to draw any conclusions about the su- 
pergravity spectrum associated to our D~particles, and in particular the topic of scaling 
solutions, in subsection 6.1. We will then discuss our most important results: the refined 
index to enumerate bound states and our computational scheme, in subsection 6.2. We 
will support our discussion with a metaphoric interpretation for the need of a refined in- 
dex for some polar states, along the lines of their interpretation as D~particle/black hole 
chromosomes given in section 2. 

6.1 Absence of scaling solutions? 

In this section, we will discuss some of our results (as well as results from [1]), and in 
particular whether they might shed any light on the mysterious nature of scaling solutions. 
We will start by reviewing some of their basic features for the reader. More details can be 
found e.g. in [3]. 

Scaling solutions for a given total charge F are microstates that owe their name to the 
special feature that the distance between their centers is not fixed, but rather a 'scaling' 
modulus. The appearance of scaling solutions can be understood easily using a concrete 
three-centered example. To make following equations transparent, the shorthand notation 
Tij = {Ti,Tj) is introduced, as well as h to label the constants appearing in harmonic 
functions, belonging to a multi-centered supergravity solution, 

^ r- 

^(^^) = E ^3-7 + ^' (6-1) 
1=1 

where h = —2lm{e~^°'^})r=o- One can find a scaling solution by treating \xi — Xj\ = XTij 
as independent variables and sending A — t- 0: this explains why these solutions are called 
'scaling solutions'. The distances between the centers are not completely independent: in 
order for such a solution to exist, one must respect the triangle inequality, 

r2i + Tis > + cyclic permutations. (6.2) 

In the limit A = (or A infinitesimally small), the locations of the centers in spacetime 
become identical, and the black hole solution becomes indistinguishable from a single- 
centered black hole for a distant observer, whereas for an observer remaining at finite 
distance from the centers, the solution stays multi-centered. The interpretation is that the 
throats of the black holes have melted together and that the near observer has disappeared 
down the throat. This is illustrated in figure 3. 
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Figure 3: Scaling solution: As the positions of several black holes come very close to each other, 
the throats belonging to the black holes melt together, eventually disappearing down the (infinitely 
long) throat of one big black hole. From the outside, it is indistinguishable from a single-centered 
black hole, for an observer placed somewhere down the throat it remains multi-centered. 

We will now discuss some microstates of the D-particles investigated in this paper 
(as well as an example from [1]), for which one would expect scaling solutions (ignoring 
regimes of validity of the description for the moment) in the supergravity approximation. 
The condition for existence of scaling solutions for a D4~D2-D0 D-particle state, (6.2), 
consisting of a fluxed D6-brane, Ti = e^(l + ^^^j^), a D6~brane, r2 = —(1 + ^^|^), and 
N DO's, Ts = (0, 0, 0, -N) reads 

• N >2 for the sextic: in particular, the charge labeled by = 0, Ago = —2 should 
allow for a scaling solution. Note that the 40'392 BPS states we found should be the 
exact number of BPS quantum states (at least it matches with the prediction from 
modularity) . 

• N > 2 for the octic: again, the charge system Aq = 0, Ago = —2 fulfills this require- 
ment. Note that the 86'140 BPS states should be exact. 

• > I for the decantic: in this case, the charge system Aq = 0,Ago = —2 should 
support scaling solutions. At this moment, we cannot address (at least predict and 
check) the number of BPS states conclusively. 

• N > ^ for the quintic, studied in [1]: this implies that the charge system Aq = 
0, Ago = —3 falls into this category. The number of 5'817'125 BPS states found in 
[1] should be exact. 

There are several possibilites: it might be, that some BPS states in the quantum 
spectrum might correspond to classical BPS states which are scaling solutions. Another 
is that scaling solutions exist, but are just not 'valid' classical microstate descriptions of 
the D-particle, as they do not correspond to any quantum states. What is meant by 
'correspond'? The authors of [17] studied the quantization of the phase space of smooth 
supergravity solutions. In well understood situations (type IIB compactifications and the 
Strominger-Vafa D1-D5-P black hole), it can be shown that the quantization of the phase 
space of smooth supergravity solutions (thus a restriction to the states carrying quantum 
numbers) is in one-to-one correspondence with the BPS microstates in the D-brane de- 
scription. Microstates on the gravity side arise as wavefunctions, which localize on a unit 
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volume of phase space. A classical solution can then be interpreted as the limit of this 
localized wavefunction. It can however be the case that there are quantum states in the 
spectrum, which do not localize on such a unit of phase space, and it is questionable if 
these states have a reasonable classical description. On the other hand, the authors also ar- 
gue, that 'nominally' classical solutions might occupy the same volume in quantized phase 
space (arise as limits of a state localized in the same volume), although they differ on a 
macroscopic scale in their gravitational description. This cannot be reasonable when taking 
into account Heisenberg's uncertainty principle, and thus 'classical solutions' of this kind 
might not be 'valid' geometries of the black hole. They might just be a peculiarity without 
physical significance, found in supergravity. This is exactly what the authors claim to be 
the case for scaling solutions. Thus, scaling solutions might not be valid classical limits of 
black hole microstates. Nevertheless, an index accounting for scaling solutions based on 
the quantization of phase space was constructed in [17], and this might help when trying 
to match classical and quantum BPS states. 

For three of the four listed charge systems, the number of BPS states in the spectrum 
can be predicted by modularity, and the methods used in this paper allow for exact results 
to be derived, based on the prescriptions for refined calculations. One is led to believe 
that one has indeed exhausted the BPS spectrum of the quantum theory, in the sense that 
one has classified and enumerated all BPS microstates of the corresponding D-particles, 
carrying the desired charge. At the same time, this means that the split attractor flow 
conjecture worked perfectly. For none of the mentioned charge systems is there seemingly 
room for states, which would be described by scaling solutions, after taking a classical 
limit. This can be concluded from the fact, that an index (as suggested in [17]) for scaling 
solutions seems to yield results that strongly deviate from any numbers contributing to the 
total number of BPS states one would expect. It seems reasonable to suppose that the exact 
numbers for BPS states of the D-particles investigated can be interpreted as an indication 
that the interpretation of [17] on scaling solutions is indeed correct. It would be wrong 
however to make an overstatement, as these implications do remain quite speculative, and 
there are no properly legitimated methods to compute an 'index for scaling solutions', 
taking into account the regime of validity needed in this case. One should at this point 
remain open to other ideas, one of which might be the wild proposal that 'special states' 
could be quantum partners of scaling solutions. The authors believe this to be rather 
unlikely, but further investigation on this exciting topic remains for the future. 

Scaling solutions might not be good classical limits of quantum microstates, or one 
might not understand them properly, yet. Maybe — being very careful not to make any 
overstatement — it is too much to hope for, that split attractor flow trees classify the full 
BPS spectrum. In any case, the split attractor flow tree conjecture is extremely accurate 
for a large number of BPS states in supergravity, and also classifies BPS bound states 
outside of the regime of supergravity: this claim was successfully put to various non-trivial 
tests in this work. 

6.2 The meiosis of D particles/black holes 

Our main results involve a scheme for computing a refined index enumerating BPS bound 
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states, and in connection with this we have calculated partitions of Donaldson-Thomas 
invariants, resulting from the distinction of constituent states as based on the perception 
of the tachyon fields gluing together a bound state. Basically, the fact that one needs to use 
a refined index to enumerate BPS bound states, results from the fact that the constituents 
are not fully independent, there is some sort of 'interaction' between the constituents for 
which we would like to give a nice little mctaphoric interpretation, using the idea to call 
polar states chromosomes of D-particles/black holes, as stated in the introduction of this 
paper. 

The degeneracy of a polar state is naively believed to factorize into three factors, 
the degeneracy of the first 'parent' constituent, the degeneracy of the second 'parent' 
constituent, and the index accounting for the tachyonic degrees of freedom, gluing the state 
together. Although this analogy again has many shortcomings, one might — when thinking 
of the suggested metaphor — count the number of possible outcomes of meiosis, which 
produces the four haploid gametes, necessary for the sexual reproduction of eukaryotic 
organisms. Naively, one might expect there to be only one outcome (with half of the 
gametes containing the father chromosome and the other half the mother chromosome). 
During the meiosis however, crossing-over (intrachromosomal recombination) takes places, 
referring to the process of exchange of pieces of the DNA on pair chromosomes (a sort 
of 'interaction' between chromosomes). This means that the number of possibilities does 
not factorize. For our analogy, the condensation of two higher-dimensional branes is the 
parallel for the reductional division of two chromosomes leading to a (haploid) cell, a 
'gamete'. For a BPS bound state, constituent branes of course do not really perform a 
crossing-over, but the number of tachyonic strings (whose specific existence does depend 
on both of the branes) gluing the state together jumps. Still, this interpretation is meant to 
stress the point which the authors believe to be important: namely, that the constituents 
of a bound state are not independent, but in general do perform some interaction. This 
can be taken into account when calculating indices by using the refined prescription, used 
in this paper. 

Finally we would like to state some implications and future directions of reseach sug- 
gested by our results: 

1. It would be extremely appealing to be able to find more certainty for the correctness 
of our new prediction of the elliptic genus for the decantic hypersurface in W^IPiii25- 
This would involve enumerating the remaining BPS states, amongst which those 
corresponding to a single flow. 

2. Probably, special loci as discussed in this paper will appear even more prominently, 
when working with Calabi-Yau varieties embedded in toric ambient spaces. As this 
feature gains more importance, it would be interesting to study tractable examples, 
and see how important this effect is. Apart from that, it would also be nice to 
predict elliptic genera for higher class divisors (also for one-modulus CY's), and also 
sec what implications the refined BPS bound state index computations have for those 
examples. A goal is to develop the prescriptions for refined bound state indices to 
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Figure 4: Polar states, non factorization and crossing over: 1. Top: the naive factorization 
of states. 2. Bottom: The number of gluing strings between two branes is not independent of the 
states in which the two branes reside. This ruins factorization of the total number of states for such 
a polar term. In meiosis, homologous pairs of chromosomes are lined up and can exchange pieces 
of their DNA (intrachromosomal recombination). This phenomenon is called crossing-over and is 
depicted for the lower pair of chromosomes. This means that the number of (possibilities to express 
the) genes, passed on to the next generation does not factorize between homologous chromosome 
pairs. 

full generality, which includes constraint loss arising from curves and particles on all 
constituents of a bound state. 

3. Additionally, it is a question what the implications are for a factorization formula for 
BPS indices, and also of the OSV conjecture. The factorization of BPS states is only 
possible on the level of densities. Thus the question is whether one can also write 
down a closed formula in the spirit of OSV for a D-particle / black hole partition 
function, on the level of index densities. 

4. Considering our interest in Donaldson-Thomas partitions, it would be interesting to 
study partition functions built from Donaldson-Thomas partitions. Given the fact 
that only three divergent asymptotic series (amongst which the Donaldson-Thomas 
partition function) approximating the topological string partition function are known, 
it would be interesting to study the convergence behavior of such a 'partition function 
of Donaldson-Thomas partitions'. 

5. Using our metaphoric interpretation of chromosomes, one might call the quest for 
finding more elliptic genera for D-particles and eventually of course black holes, the 
black hole genome project. It would seem an interesting task to map mixed ensembles 
to IIB string theory, lift them to five dimensions, and in this way strengthen the link 
between the formalism of 4d polar states (the chromosomes of black holes), and the 
fuzzball program in 5d. A first step in the direction of mapping polar states to 
fuzzball microstates in 5d was taken in [29], for BPS polar states in A/" = 8 (4d) 
supergravity theories. 
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A. Setup and conventions 

As all Calabi-Yau 3-folds X used in this paper are one-modulus examples, h^^'^\X) = 1, 
we use H to denote the basis element of H'^{X,Z). We also use the symbol H^*{X,Z), or 
for short, H'^*{X) to denote the even degree cohomology of X: 

H^*(x) = H°{x) e H^{x) e H'^ix) e h^{x). (a.i) 

For a type IIA string theory compactification, a vector bundle with fibers H^* (X) over the 
Kahler moduli space is defined, by virtue of the properties of special geometry. These fibres 
are naturally endowed with a symplectic product (Dirac-Schwinger-Zwanziger) , (•,•) : 

(ri,r2) = / FiAF^ (A.2) 

Jx 

where it is understood that only the wedge products of the right dimension remain as an 
integrand, and r2 is defined under an operation inverting the sign of the two-form and the 
six-form components: F = pO + T^H + T'^H^ + T^H^ ^ r* = F^ - T'^H + T^H^ - T^H^. 

The IIA complexified Kahler form will be written as^ tH = {B + iJ)H using our basis, 
and allows us to form the holomorphic period vector 

Oioi := -eMtH) = -l-tH- ^i?^ _ ^^3 

The Kahler potential of the Kahler moduli space reads 

K = -In (i(J^hol,nhoi)) = -In Q J^iJ A J A J) 

which also allows us to define the normalized period vector 

:=exp(^)J7hoi- (A.3) 

The general formula for the induced charges on a D-brane (due to the WZ term in the 
Born-Infeld action) wrapped on a (sub)-manifold W reads 



cDbrane n — 



Jw V A{NW) 



^We will also refer to t just as the Kahler modulus. 
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where A is the A-roof characteristic class, TW the tangent bundle of the brane, and NW 
its normal bundle. From this formula, one obtains for a D6~brane carrying U(l)-flux with 
field-strength Fi, a D2 of class — /3i and A^i DO's the following polyform: 

rz)6 = (l - /3i - x{Cp,) + iVi) (l + ^) , (A.5) 

where /3 G H^{X^'L), and C2(X) is the second Chern class of the tangent bundle of the CY 
threefold X. Similarly, a D6 with flux F2 will bind to a D2 of class /32 and DO's to give 
the following total charge vector: 

= (1 - /32 + {\ ) + N,) c) (1 + ^) , (A.6) 

The modification with respect to the general formula is the addition of D2 and DO charge 
in the form of sheaves, which can be thought of as generalizations of bundles (U(l) fluxes). 
Notice that the D6 will bind with a D2, the D6 with a D2, but both will bind to DO's. 

For a D6-D4-D2-D0 brane system, we denote the charges as {pQ,p,q,qo), such that a 
polyform using our basis can be written as 

T=po+pH+^H^ + ^H^ (A.7) 

where H := fj^ H'^. We will equivalently use a vector notation: F = {po,p, q, go)- 

The central charge of a D-brane system F is defined using the period vector (A. 3), 



Z = (F, n) = exp(-) ^po-f' - p-t' +qt- qo j , (A.8) 

The central charges of brane systems under investigation will be used for establishing the 
existence of single/split flows and thus of BPS states, as will be explained in subsection C. 

B. Modified elliptic genera and D particles 

Modified elliptic genera arise as supersymmetric BPS indices of the (0, 4) Maldacena- 
Strominger-Witten CFT [30]. A non-linear sigma model realization of this CFT has been 
derived from the dimensional reduction of the M5-brane worldvolume perspective. In this 
paper, we will work from the II A perspective. In the following, we use a basis G 
H4^{X,'Ij) and define the intersection numbers as QDabc- The definition Dab = Dabcp'^ 1 
its inverse and D = DabcP^ 

p^p'-' will also come in handy. Using mixed ensembles 
of D4-D2-D0 branes with fixed magnetic D4-brane charge p^ and variable electric charges 
(qa, Qo), the central charges of the CFT read 

CL = P^+C2-P, Cjf = P3 + • P, (B.l) 

where C2 is the second Chern class of the Calabi-Yau X and P is the Poincare dual of the 
hypersurface p. The modified elliptic genus is defined as 

Z{q,q,y) = Ttr ^^^(-1) V"" - 9"^''- - e^"^^"^ j , (B.2) 
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and is constrained to be a weak Jacobi form of weight (— |, |)- The flux on the D4-brane 
is an element of H^{P,'Z) + ^, where the ^ factor accounts for the Freed-Witten anomaly. 
Each flux can be decomposed as follows 

i^ = ^ + /"+/^+7, (B.3) 

where /H E Lx = i*H'^{X,7j) and /"*" G Lj^, and 7 is given by gluing vectors in the flux 
lattice. The reader can consult [3] for more details. The modifled elliptic genus can be 
decomposed as 

Z{q,q,y) = ^Z^{q)e^{q,q,y) , (B.4) 

7 

where 7 runs through a finite set of gluing vectors. The Z^{q) are meromorphic functions 
of the variable q and the theta functions are so-called Siegel-Narain theta functions: 

e,{q,q,y) = E (-ir"''-g^('^'''"')'"'"r^^^"^^^'e2-^"'^- (B.5) 

where the summation over qA is understood to be qA = 6Dab{^P^ + k^) + 7A, & 
and n is the dimension of the flux lattice. 

The dimensionality of the vector Z^ is thus given by the number of independent el- 
ements 7 of the discriminant group, the gluing vectors. Additionally, as discussed in 
[4, 5, 3, 31], it follows from modular invariance, that there is an identification 

Z-y = Zs, for 7 = — (5 mod Lx- (B.6) 

The Z-y's can be given some physical interpretation: Zq corresponds to a sum over 
states with no added D2-charge, and increasing DO-charge as the powers of q increase. 
Each coefficient in the g-expansion corresponds to the index of a state with fixed DO- 
charge. Similarly, the Z^'s correspond to states with added D2-charge^. Schematically, 
the first few terms of Z^ will look as follows 

^7(9) = + # + # 9' + •••) , (B.7) 

where % = qQ — j^D^^^qAqB^ and a is the highest possible value of ^0 for a given 7. In 
this paper, only the Z-y{q) functions will be of concern, as they contain all the relevant 
information. A stringent mathematical property of weak Jacobi forms is the fact that they 
are entirely determined by their polar part, i.e. terms with negative powers of q. These 
terms correspond to charge configurations that satisfy % = qo — ^2^"^^^^^^ > 0. Such 
configurations will be referred to as polar states. To construct elliptic genera, we used the 
method of generating modular representations as in the appendix of [2] . 

®Note that by spectral flow, which transforms qA — >■ <1A — BDAsk^ and go — (Jo + k^QA — SDAsk^k^ , 
most of the D2-charge can be absorped in qo 
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C. The split attractor flow tree conjecture 



In the following, the reader will find a very brief presentation on the use and meaning 
of split flow trees. By the split attractor flow tree conjecture, [3], single flows and split 
flow trees are believed to be an existence criterion and provide a complete classification^ 
for BPS states in type II string theory. Single/split flows are graphical depictions of the 
flow of the Kahler moduli (or complex structure moduli in the mirror type IIB picture) 
belonging to a BPS solution of supergravity, but their meaning extends beyond the range 
of validity of the supergravity approximation to type II string theory, as suggested by the 
smooth interpolation between the supergravity picture and the quiver description of D- 
brane bound states found in [6]. 

In order to incorporate the use of split attractor flow trees also outside of the regime 
of validity of the supergravity approximation to string theory, the following definition of a 
split flow tree is given: 

A split flow tree belonging to the total charge T and n constituents (ri,...,r„) is a 
set of data 



consisting of 2n points, a background value too, a set of ?i — 1 split points tj,spiit (j = 
{l,..,n — 1)), and n attractor points tj* (j = l,...,n), one for each center. Using this 
notation, a single flow, for example, is denoted as (too;^*)- Note that this definition boils 
down to the most essential features, but the intuition, that the flow tree depicts the values 
of the scalars belonging to a supergravity solution in moduli space, remains valuable, when 
constructing split flow trees. 

The construction of a split flow tree 

A split flow tree is built as follows. One follows the incoming branch^ of a flow tree 
from radial infinity towards a putative attractor point, until one hits a wall of marginal 
stability for two non-local^ constituents (ri,r2) / such that F = Fi + This is de- 
fined as the hypersurface in moduli space where the phases of the central charges align, 
arg(Zi) = arg(Z2). The modulus of a central charge measures the mass of the correspond- 
ing state, whereas the phase indicates which J\f = 1 supersymmetry of the original M = 2 
supersymmetry is preserved by the state. If the phases of two central charges align, the two 
states are mutually BPS (i.e. preserve the same supersymmetry) and the binding energy 

^The statement that they classify all BPS states in string theory is sometimes referred to as the strong 
form of the conjecture, while one can also formulate the weaker conjecture for BPS states in supergravity. 

*The term 'incoming branch' refers to the part of the flow tree connecting the background point and the 
first split point. 

®If one starts with local constituents, (ri,r2) = (such as two branes of the same type), they are 
trivially mutually BPS, and one speaks of a threshold wall instead of a wall of marginal stability. 
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of the BPS bound state vanishes, 

IZ1+2I = |Zi| + IZ2I, (C.2) 

or equivalently, 

Re(ZiZ2) > 0, Im(ZiZ2) = 0. (C.3) 

If one reaches the waU from the side where (Fi, r2)(arg(Zi) — arg(Z2)) > 0, the decay of 
r — >■ Fi + r2 is energeticahy favored. MicroscopicaUy, the non-zero value of (Fi, r2) means 
that there are chiral strings stretched between the two constituent branes, which make the 
decay (or the recombination if one goes in reverse) possible. One then follows the flows of 
the constituents, which might decay again according to the same scheme, until every end 
branch flows towards an attractor point. This is illustrated by the following figure for a 
three-centered solution: 



Figure 5: Split flow tree for a three centered solution with charges F = Fi + V2a + ^2b'- 
One starts at the background point, too, at the top, and follows the incoming branch (plotted in 
black) until one hits a wall of marginal stability (plotted in red) between Fi and F2. This is the first 
split point t^g"^ , wherefrom single flows corresponding to the two centers continue (both branches 
are plotted in red). The first center reaches an attractor point, il^-*, whereas the second branch 
again reaches a wall of marginal stability (plotted in green) at t^g'"^^. At this point, the charge F2 
splits into F2a and F2f,. Two branches, corresponding to these two charges, are plotted in green: 
they flow to their attractor points, tl ' and t\ ' . 

For a given total charge F one might in general find several different split flow trees, 
and maybe also a single flow, all contributing to the total index of BPS states with this 
total charge. In general, the BPS spectrum of states ri(F) remains invariant under in- 
finitesimal variations of the background moduli too, but it can jump when the moduli are 
driven through a wall of marginal stability. This is intuitively clear; a certain split might 
either not be possible anymore because one can only reach the appropriate wall of marginal 
stability from the unstable side, or, alternatively, a new type of split becomes possible, as 
one can now reach this wall from the stable side. In such a case one has taken the back- 
ground into a different area, and the fact that there are different basins of attraction in 
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moduli space has led to the name area code for the backgromid. Hence, an index of BPS 
states should more precisely be denoted as f](r,too)- 



Threshold walls 

Apart from the walls of marginal stability, which separate regions in moduli space be- 
tween which a BPS index can jump, there is a second type of wall that is of importance 
for our work: walls of threshold stability. The distinction between these two kinds of walls 
is explained in more detail in [15]. The threshold conditions are the same as for marginal 
stability 



however this time, the two charges are mutually local (ri,r2) = 0. Microscopically this 
means that there are no tachyonic strings between the two branes that can condense to 
merge Fi with Most importantly, a BPS index cannot jump when crossing the thresh- 
old stability wall with the background modulus. What can happen is that the flow tree 
changes topology: one can imagine, for example, a three-centered solution, with one 'satel- 
lite' center bound to either of the two others, which then changes sides and is bound to 
the other center after crossing the threshold wall with the background modulus. 

This allows the formulation of the split attractor flow tree conjecture ([3]): 

The split attractor flow tree conjecture (strong version) 

1. For a given background too in Kahler/complex structure moduli space, the existence 
of a split attractor flow tree starting at the background too, with a given total charge 
F and endpoints corresponding to Fj, is equivalent to the existence of a BPS bound 
state in string theory, with constituents Fj. 

2. The number of split flow trees and hence the total number of states with a given charge 
F in a fixed background is finite, at least when charge quantization is imposed. 

As will be discussed later on, our results offer support for a strong version of the conjecture. 

In this work, we use split flow trees to establish existence of BPS states in a regime where 
supergravity cannot be trusted. As the D-brane charges used are very low, the attractor 
equations will drive the horizon size to be very small with respect to the string scale. Thus, 
the supergravity approximation breaks down, and higher curvature corrections from per- 
turbative a' effects become important. At the same time, the cycles on which the branes 
are wrapped also become very small, and worldsheet corrections to the central charges of 
the considered brane systems become dominant. Luckily, one can calculate exact central 
charges in the mirror picture, as the periods are exactly determined by classical geometry. 
The periods are determined as solutions to the Picard-Fuchs differential equations and 
cannot be written down and used analytically. Therefore, we used numerical techniques to 
tackle this problem. These techniques were developed in [10], and are partially based on 





(C.4) 
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the techniques prepared in [32]. They were also used in the same way as in this work in 
[1]. The reader will find more explanations in appendix D. 

D. Methods to establish split flows outside the supergravity regime 

For configurations with low D-brane charges, the attractor flow equations drive the horizon 
size of solutions to very small sizes in string scale units. This automatically leads far outside 
the supergravity regime, requiring higher curvature corrections. However, as stated in the 
main text, the main tool of analysis, namely split attractor flow tree techniques, retains its 
meaningfulness as predicted by the strong split attractor flow conjecture. As the attractor 
equations will typically drive the cycles and the CY itself to stringy sizes, the central charges 
of brane systems receive important worldsheet instanton corrections as a consequence, and 
these need to be taken into account. In type IIA string theory this would be an impossible 
task. Luckily, mirror symmetry solves the problem as the central charges are exactly 
determined classically by the periods of the holomorphic three-form in type IIB string 
theory on the mirror CY manifold. The mass of a BPS saturated brane wrapping an even 
cycle corresponds exactly to the notion of quantum volume of that cycle from [32]. The 
general scheme is thus to identify an integral basis of three-cycles, calculate the periods of 
these cycles, find the explicit mirror map and in this way define the quantum volume of 
any even dimensional cycle of the mirror. 

The periods of the mirror Calabi-Yau's are determined as solutions of the Picard-Fuchs 
equation, and are known as Meijer G-functions. Thus, it is not possible to write down 
analytic formulae for central charges, and certainly not possible to work with split attractor 
flow trees analytically. This hurdle is overcome in the following by means of numerical 
approximations with Mathematical*^. The periods of the mirror CY are evaluated, and a 
lattice of points is created, from which the function can be approximated by interpolation. 
Split flow trees and single flows are then established numerically. More details on the 
applied technique can be found in [10]. The mirror symmetry induced monomial-divisor 
map is used to convert Kahler to complex structure modulus, and map {D6, D4, D2, DO) 
brane systems Fa into their {D3, D3, D3, D3) brane mirrors F^, L : H'^*{X,Z) H-^{Y,Z), 
and then analyze the attractor flows of the exact central charges in complex structure 
moduli space (or more precisely in the n-fold cover w-plane). In the IIB picture, BPS 
black holes are made of D3-branes wrapped along special Lagrangian three-cycles of the 
internal Calabi-Yau manifold, whereby these D3's (and their corresponding three-cycles) 
can split up into intersecting D3's, by moving in the complex structure moduli space of the 
CY across some 'line of marginal stability'. 

The mirror map L can be found from a comparison of the IIA and IIB periods near the 
Large Complex Structure (LCS) point. It relates the even cycles of real dimension 2j on 
the CY Xa to the three-cycles 7j of the mirror CY manifold Xb- The periods Ilj = 0, 
of the holomorphic three-form on the 7^ have leading logarithmic behavior log-'(z) near 
z = (LCS point), using the coordinate z = ip~^ and ■ip~^^ on the mirror moduli 

^"Special thanks to F. Denef for sharing his Mathematica code, [10]. 
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space, for the sextic, the octic, and the decantic, respectively. The periods are solutions to 
the Picard-Fuchs equation, 



zd, H {zd, + Pi-l)-z H {zd, + aj)\u = 0, (D.l) 

i=l...q j=i-P 

where the aj and the /3j are model dependent constants, which read /3j = 1, i = 1,2,3, 
and a = (g, ^, |, |) for the sextic, (|, |, |, |) for the octic, and {j^, j^, ^) for the 
decantic hypersurface, used in this paper. The Meijer-functions (G-functions) Uj{z) can 
be expressed as 

This particular basis of periods is related to three branching points (LCS point, conifold 
point and Gepner point) which are connected by appropriately chosen branch cuts. 

Associated to these branch cuts are three types of monodromies, which can be ex- 
pressed as matrices acting on the periods. The monodromy T(0) around the LCS point 
(z = 0, = 00) and the monodromy T'(oo) around the Gepner point {z = oo, ^ = 0) act 
on the period vector U(z) = {Uj{z))j=i,_A as follows, 

\]{e^^^z) = r(0)U(z), \z\ « 1 

U(e2'^'z) = r(oo)U(z), \z\ » 1. (D.3) 

The third monodromy matrix of course follows directly from the other two, as a monodromy 
can always be seen either as 'around one of the branching points' or, equivalently as a 
monodromy 'around the two other branching points' in the appropriate directions. For the 
monodromy around the conifold point one has 

r(l) = r(cx))T(0)-^ lm{z) < 0, 

r(l) = r(0)-^T(oo) lm{z) > 0. (D.4) 

Using Mathematica, we evaluated the periods for our models on a lattice of points, the 
number of which can be adjusted to the degree of precision demanded. By interpolation 
this gives a numerical approximation of the periods. 
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